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Abstract 

We present a QCD sum-rule determination of the heavy-quark kinetic 
energy inside a heavy meson, —\\/2rriQ, which is consistent with the field- 
theory analog of the virial theorem. We obtain — Ai ~ (0.10 db 0.05) GeV 2 , 
significantly smaller than a previous sum-rule result, but in good agreement 
with recent determinations from the analysis of inclusive decays. We also 
present a new determination of the chromo-magnetic interaction, yielding 
\ 2 (m b ) = (0.15±0.03) GeV 2 . This implies m|* -m% = (0.60±0.12) GeV 2 , 
in good agreement with experiment. As a by-product of our analysis, 
we derive the QCD sum rules for the three form factors describing the 
meson matrix element of a velocity-changing current operator containing 
the gluon field-strength tensor. 
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1 Introduction 



The physics of hadrons containing a heavy quark simplifies greatly in the limit 
where the heavy-quark mass mg is taken to infinity. Then new symmetries of 
the strong interactions arise, which relate the long-distance properties of many 
observables to a small number of reduced hadronic matrix elements 1|. A 
systematic expansion around the heavy-quark limit has been applied success- 
fully to learn about the properties of heavy mesons and baryons, such as their 
spectroscopy and decays. 

A convenient tool to study the implications of the heavy-quark limit and to 
perform the 1/mq expansion is provided by the Heavy-Quark Effective Theory 
(HQET) 0, which is constructed by introducing a velocity-dependent field h v (x) 
related to the original heavy-quark field Q(x) by 

K(x) = e m ^ x ^Q(x), (1) 

so that j)h v = h v . Here v is the 4- velocity of the hadron containing the heavy 
quark. The phase redefinition in (|]) removes the large "mechanical" part rrigv 
of the heavy-quark momentum, which is due to the motion of the heavy hadron. 
The field h v (x) carries the "residual momentum" k = pq — rriQV, which arises 
from the predominantly soft interactions of the heavy quark with gluons. The 
effective Lagrangian of the HQET is [|-[[K| 

£ eff = h v iv-D h v + —^—h v (z-D±) 2 ^+C ma g(m Q //i) h v a„ v G^h v +0(l/m 2 Q ) , 
2m Q 6 4m Q ^ v 

(2) 

where = d^ — ig s A^ is the gauge-covariant derivative, and Dj_ = D^ — (v-D) v^ 1 
contains its components orthogonal to the velocity. The gluon field-strength 
tensor is defined as [D^,D U ] = —ig s G^ v . The origin of the operators arising at 
order 1/m.Q in (0) is most transparent in the rest frame of the heavy hadron: 
the first operator corresponds to the kinetic energy resulting from the residual 
motion of the heavy quark inside the hadron (note that (iD±) 2 = — (2D) 2 in the 
rest frame), whereas the second operator describes the magnetic interaction of the 
heavy-quark spin with the gluon field. The Wilson coefficient C mag results from 
short- distance effects and depends logarithmically on the heavy-quark mass and 
on the subtraction scale /x, at which the chromo-magnetic operator is renormalized 
0. As a consequence of a reparametrization invariance of the HQET, the kinetic 
operator is not renormalized \TT\. 



In many phenomenological applications of the HQET, the forward matrix el- 
ements of the dimension-5 operators in (0) play a most significant role. They ap- 
pear, for instance, in the spectroscopy of heavy hadrons [fO^— IP-4|, in the descrip- 



tion of inclusive decay rates and spectra |15||-|pl]|) as well as in the normalization 



of transition form factors at zero recoil [IB]. For the ground-state pseudoscalar 
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and vector mesons, M = P and V, one defines two hadronic parameters, Ai and 
A 2 (/i), byfl 

(M(v)\h v {iD x fh v \M{v)) = Ai, 
(M(v)\ h v g s o^h v \M(v)) = 2d M \M , (3) 

where we use a mass- independent normalization of states such that (M(v) \ h v h v \M(v)) 
1. The coefficient d M takes the values dp = 3 and dy — —\ for pseudoscalar 
and vector mesons, respectively. The \x dependence of the parameter A2 cancels 
against the \i dependence of the coefficient C mag in (fj). The product C mag A2 
is renormalization-group invariant. This quantity can be extracted from spec- 



troscopy using the relation [13 

i (m 2 B , - ml) = C mag (l) X 2 (m b ) + 0(A 3 /m b ) = 0.12 GeV 2 . (4) 
In leading logarithmic approximation, one finds |9j 

CWmcj/M) = ( ^^V //3 ° , (5) 



a s (/z) 



where /?o = 11 — is the first coefficient of the (3 function. Using this result, 
we obtain 

\ 2 {m b ) ~ 0.12 GeV 2 , A 2 (/i ) ^ 0.15 GeV 2 , (6) 

where we have used the values a s (m&) = 0.21 and a s (/xo) = 0.4, corresponding to 
a low renormalization point /^o ~ 1 GeV. 

Although spectroscopic relations may be used to extract the differences of 
the matrix elements of the kinetic operator between different hadron states, the 
value of the parameter Ai itself cannot be determined from spectroscopy. Indeed, 
at present it is not even known whether Ai is a "physical" parameter in the 
sense that it can be defined unambiguously in a non-perturbative way. It may 
be that any such definition is intrinsically ambiguous because of the presence 
of infrared renormalons; however, it may also happen that for some fortuitous 



reason renormalons are absent in the case of the kinetic operator |22|, |23|. From 
the practitioner's point of view, Ai becomes a useful parameter once a scheme for 
treating perturbative corrections in the HQET has been specified. Here, we shall 
work in the Ms subtraction scheme and assume that our value of Ai is used in 
connection with theoretical expressions that have one-loop perturbative accuracy. 

Theoretical information about the parameter Ai can be obtained from a non- 
perturbative evaluation of the first matrix element in (|3|). Using QCD sum rules 
in the HQET, Ball and Braun have derived the value -\ x = (0.52 ± 0.12) GeV 
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24}| , which has been adopted subsequently in many phenomenological analyses. 



1 Another common notation is to define /i^ = — Ai and /jIq = 3X2- 
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This value is surprisingly large; it implies an average momentum of the heavy 
quark inside the meson of order 600-800 MeV. In fact, an earlier QCD sum- 
rule calculation using a less sophisticated approach had given the smaller value 
— Ai = (0.18 ± 0.06) GeV 2 |j25| . A theoretical argument in favour of a smaller 
value of the kinetic energy was presented in Ref. [Q: the field-theory analog of 
the virial theorem relates the first matrix element in (|3|) to a matrix element of an 
operator containing the gluon field-strength tensor, making explicit an "intrin- 
sic smallness" of Aj. On the other hand, a large value of — Ai was argued for by 
Voloshin and by Bigi et al., who derived the lower bound — Ai > 3A2 — 0.36 GeV 2 
using first a quantum-mechanical reasoning J27j and later field-theoretical argu- 
ments based on zero- recoil sum rules ||28|| . Recently, however, it was shown that 



this bound is weakened significantly by higher-order perturbative corrections [^9 
Thus, small values of — Ai can no longer be excluded a priori. Several authors 
have attempted to extract Ai (together with the "binding energy" A) from a 
combined analysis of inclusive decay rates and moments of the decay spectra in 
beauty and charm decays |3(J-[[£|. The most recent of such analyses give val- 
ues —Ax w 0.1-0.2 GeV 2 . A value of Ai has also been extracted using a lattice 
simulation of the HQET, yielding -X x = (-0.09 ± 0.14) GeV 2 Q. 

In view of these developments, it seems worth while to reconsider the problem 
of calculating Ai and A2 using QCD sum rules. Here we present the results of 



a new analysis, which is based on the virial theorem JIB], |26|]{] Our approach 
has the advantage that both parameters are determined simultaneously from the 
zero-recoil (i.e. equal- velocity) limit of the QCD sum rule for the matrix element 
of a local dimension-5 operator between two meson states moving at different 
velocities. As a result, Ai and A 2 are obtained from a set of two sum rules with 
very similar systematic uncertainties. 



2 Derivation of the Sum Rules 

The central object of our study is the meson matrix element of a local dimension- 
5 operator containing two heavy-quark fields at different velocities. Using the 
covariant tensor formalism of the HQET |36 |, we write 

(M'{v')\ h v ,Tig s G^h v \M(v)) = -Ti {^(v,v')M' (v')TM(v)} , (7) 

where T is an arbitrary Dirac matrix, and 

, j, x 1 + v f 75 ; pseudoscalar meson P(v) 
M(v) = — -r { j . i/v \ ( 8 ) 

2y 2 I f ; vector meson V (v, e) 

is a matrix representing the spin wave- function of a ground-state meson M moving 
at velocity v. The most general decomposition of the tensor form factor <ft iV 



Preliminary results of this analysis have been presented in Ref. 35 
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consistent with Lorentz covariance and heavy-quark symmetry reads 



4T{v,v') = 0i H {vN' v -v v v'^)+(j) 2 {w) [(v - vYY -(v- v'Yl^+hH ^ , 

(9) 

where w = v ■ v' . For equal velocities (w = 1) only the last term appears, and 
comparing with (H) we obtain the normalization condition 

3 (1) = A 2 . (10) 

The equations of motion of the HQET imply another normalization condition for 



a certain combination of the invariant functions 4>i{w). It reads [13 



30 1 (1)-30 2 (1)-^ 3 (1) = -A 1 . (11) 

This relation is remarkable in that it relates the matrix element of the kinetic 
operator in ([3|) to a matrix element of the gluon field-strength tensor, in accor- 
dance with the picture that the residual motion of the heavy quark inside the 
meson is caused by its interactions with gluons. Eq. (O) can be interpreted as 
the field-theory analog of the virial theorem, which relates the kinetic energy to 
a matrix element of the "electric" components of the gluon field . 

We shall now derive the Laplace sum rules for the invariant functions 4>i(w). 
The analysis proceeds in complete analogy to that of the Isgur-Wise function. 
For a detailed discussion of the procedure and notations, the reader is referred to 
Ref. [07[| . We consider, in the HQET, the 3-point correlation function of the local 



operator appearing in (|7|) with two interpolating currents for the ground-state 
heavy mesons: 

dxdye^ k '- y - k x) ( | T {qT M ,h v ,{y), h v ,T ig s G^h v (0), h v T M q(x)} | ) 

= Tr{<^>y,MOrM'^r^r M j , (12) 

where k and k' are the external momenta. Depending on the choice = 75 or 
Tm = 7° — v a , the heavy- light currents interpolate pseudoscalar or vector mesons, 
respectively. The Dirac structure of the correlator, as shown in the second line, 
is a consequence of the Feynman rules of the HQET. The quantity obeys a 
decomposition analogous to @, with coefficient functions $j(u;, u', w) that are 
analytic in the "residual energies" uj = 2v ■ k and uj' = 2v' ■ k', with discontinuities 
for positive values of these variables. These functions also depend on the velocity 
transfer w = v ■ v ' . 

The idea of QCD sum rules is to relate a theoretical approximation to the 
Operator Product Expansion (OPE) of the above correlator to a hadronic rep- 
resentation of the same correlator in terms of physical intermediate states. The 
lowest-lying states are the ground-state mesons M(v) and M'(v') associated with 
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the heavy-light currents. They lead to a double pole located at u = uj' = 2A, 
where A = uim — rriQ is the "effective mass" of the ground-state mesons in the 
HQET [|^|. The residue of this double pole is proportional to the invariant 
functions <fii(w). We find 

^poic, , s A^H F 2 

^ (w ' w ' w) = F^TW^j' (13) 

where F is the meson decay constant in the HQET (F ~ fu^^u)- In the 
deep Euclidean region the correlator can be calculated perturbatively because 
of asymptotic freedom. The main assumption behind QCD sum rules is that, 
at the transition from the perturbative to the non-perturbative regime, confine- 
ment effects can be described by including the leading power corrections in the 
OPE. They are proportional to vacuum expectation values of local quark-gluon 



operators, the so-called condensates Following the standard procedure, we 



write the theoretical expressions for $j as double dispersion integrals and per- 
form a Borel transformation in the variables uj and uj'. This eliminates possible 
subtraction polynomials and yields an exponential damping factor in the disper- 
sion integrals, which suppresses the contributions from excited states. Because 
of heavy-quark symmetry, it is natural to set the associated Borel parameters 
equal: r — r' = 2T. Following Refs. [|37| |4(|, we then introduce new variables 
u + = \{oj + uj') and u;_ = uj — uj' , perform the integral over uj-, and employ 
quark-hadron duality to equate the remaining integral over u + up to a "contin- 
uum threshold" uj c to the Borel transform of the double-pole contribution in fllS]) . 
This yields the Laplace sum rules 

<j>i(w) F 2 e- 2K ' T = J du + e~"+' T p t (u + , w) . (14) 
o 

The spectral densities pi(u + ,w) arise after integration of the double spectral 
densities over c<j_. 

As pointed out above, the theoretical expressions on the right-hand side of 
the sum rules consist of perturbative and non-perturbative contributions. The 
leading terms in the OPE arise from the diagrams shown in Fig. |T|. In our anal- 
ysis, we shall include the non-perturbative contributions of the quark conden- 
sate (qq), the gluon condensate (a s G 2 ), and the mixed quark-gluon condensate 
(q g s o~nvG^ v q) = ml (qq). For a consistent calculation at order a s , we calculate 
the Wilson coefficients of the quark and gluon condensates to one-loop order, 
and the coefficient of the mixed condensate at tree level. At higher orders in 
the OPE, one encounters a proliferation of condensates whose values are essen- 
tially unknown. The terms of dimension six, in particular, consist of four-quark 
and three-gluon condensates. For an estimate of such contributions we include 
the effects of four-quark condensates, which arise from the diagram shown in 
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to 





to 



Figure 1: Non- vanishing diagrams for the 3-point correlator: (a) perturbative 
contribution, (b) quark-condensate, (c) gluon-condensate, (d) mixed-condensate, 
and (e) four-quark condensate contributions. The velocity-changing current op- 
erator is denoted by a white square, the interpolating meson currents by gray 
circles. Heavy-quark propagators are drawn as double lines. 



Fig. |T](e) . The calculation of the condensate terms is most conveniently per- 
formed using the fixed-point gauge x ■ A(x) = with the origin chosen at the 
position of the velocity-changing heavy-quark current. The most complicated 
part of the calculation is, however, to evaluate the perturbative contribution of 
the two- loop diagram shown in Fig. [l|(a). We have calculated this diagram using 
the techniques developed in Ref. ET|. Our results for the Laplace sum rules are: 



^(w)F 2 e- 2A / T 
^(w)F 2 e- 2K ' T 



2a„T* 



7T 3 


\W + lJ 


2a s T 


5 2 




w + 1 


T 




487T 


w + 1 ^ 


4a,T 5 





5 4 (u c /T) , 
5 A {u c /T) + (qq) 5 1 (u c /T) 



3tt 



(a s G 2 ) 5 (u c /T) + 



12T 



+ 



7T 3 W + 1 

T 



5 4 (uj c /T) - 



8a s T 2 
3tt 



(a s G 2 ) S (cu c /T) - 



(qq) 5i(u c /T) 
m 2 (qq) (0 6 ) 



(15) 



24tt w + 1 " a ~ ' ~ vv " w ' 12 12T 
The functions S n (u c /T) arise from the continuum subtraction and are given by 



1 x 

5 n (x) = — dtt 
n\ J 



n e -t 



-x ST^ 

k=0 



(16) 
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The four-quark condensate (Qq) is defined as 



— 1 f)7T 

(Qe) = 9 2 s (<tft a q E iW> = — 5- « «s(gg) 2 • (17) 
/ J 

Assuming factorization of the four-quark operator [|39| corresponds to setting 
K = 1. 

In the next step, we evaluate the sum rules in (p~5|) for u> = 1 and use the 
normalization conditions (TlO| ) and ( [TT]) to obtain the Laplace sum rules for the 
hadronic parameters Ai and A2. This leads to 



A, F 2 e~ 2 ^ = ^ 5^JT) + <M (1 - i 



A 2 F 2 e~ 2K ' T = ^ 5,{u c /T) - (qq) 5 1 (u c /T) 

7T 07T 



where 



+ ^ (* S G 2 ) 5^ C /T) - (1 + £6 ) , (18) 

e 6 = -^.^M. (19) 

For all reasonable values of the parameters, Eq is of order a few per cent, which is 
much less than the uncertainty in the parameter rriQ. Therefore, the contribution 
of the four-quark condensate can be safely neglected in the numerical analysis, 
and we shall set Eq = hereafter. The sum rule for A2 (without the contribution 
of the four-quark condensate) coincides with the result derived by Ball and Braun 
H?3|j ; our sum rule for Ai is new. Notice that the sum rule for Ai does not receive 
contributions from the quark and gluon condensates. This is a consequence of the 
fact that (in the fixed-point gauge) the light quark interacts only with the mag- 
netic components of the gluon field [03 , whereas the combination of form factors 



defining Ai in (|TT|) corresponds to a matrix element of the electric components 
0- 

For the evaluation of the sum rules, it is convenient to eliminate the explicit 
dependence on the parameters F and A by using the well-known sum rule for the 
correlator of two heavy-light currents |43|, ^37], fPJ :[] 

F i e -*,T = ^ S2MT) _ m+ Tm. (20) 



Dividing the sum rules in fll8|) by the sum rule in (|20|), we obtain expressions 
for the parameters Ai and A2 as functions of the Borel parameter T and the 
continuum threshold u c . This procedure reduces the systematic uncertainties in 



3 Since the leading terms in ( |l5| ) are proportional to a s and we have no control over the 
0(a 2 s ) corrections, we do not include the known 0(a s ) corrections in (|20|) for consistency. 
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the calculation. Moreover, it eliminates the dependence on the parameter A, 
which is known to suffer from a renormalon ambiguity problem [f22 , 45] . For 



the QCD parameters entering the theoretical expressions, we take the standard 
values |39j 

(qq) = -(0.23±0.02) 3 GeV 3 , 
(a s GG) = (0.04 ±0.02) GeV 4 , 

m 2 = (0.8 ± 0.2) GeV 2 , (21) 

as well as a s = 0.4. These values refer to a renormalization scale /xq — 2A « 1 
GeV, which is appropriate for evaluating QCD sum rules in the HQET. We shall 
comment below on the sensitivity of our result to the choice of the condensate 
parameters. 

The sum-rule parameters u c and T should, in principle, be determined in a 
self-consistent way by requiring optimal stability of the results under variations 
of the Borel parameter inside the region where the theoretical calculations are 
reliable. For too small values of T, the OPE diverges, whereas for large values of T 
the contributions to the sum rules from higher resonance states become more and 
more important. Unfortunately, the continuum-contamination problem is rather 
severe in the case of sum rules for the matrix elements of higher-dimensional 
operators such as Ai and A 2 . This is exemplified by the leading perturbative 
contribution to the correlator in (|12"D, which is proportional to 



4!T 5 



J doo + u 4 + e" w+/T = 5 4 (uj c /T) + [l - 5 A {u c /T)] . (22) 
o 



The first term on the right-hand side is assigned to the ground-state, whereas 
the second term is removed in the continuum subtraction. For the central values 
T = 0.9 GeV and uj c = 2.0 GeV determined below, the ground-state contribution 
is only 7.5% of the total perturbative contribution. For comparison, in the case 
of the sum rule in (p0|), the leading perturbative contribution to the correlator 
is proportional to 52(w c /T) + [1 — <5 2 (^ C /T)], and the ground-state contribution 
amounts to 38%. For this reason, it is better to determine the allowed regions for 
the parameters u c and T by requiring stability of the sum rule (|20| ) for the meson 
decay constant, and then to use the same values in the evaluation of the sum 
rules for Ai and A 2 [ 24fl . One finds that uo c = (2.0 ± 0.3) GeV, and the "stability 



window" for the Borel parameter is 0.6 GeV < T < 1.2 GeV [37|, f|4] . 

Using these ranges of parameters, together with the central values of the 
condensates given in (pT|), we obtain the results shown in Fig. ^|. The sum rule 
for the parameter describing the chromo-magnetic interaction of the heavy quark 
exhibits very good stability. Taking an average over the sum-rule window, we 
obtain 

A 2 (/i ) = (0.19 ± 0.02 ± 0.02) GeV 2 , (23) 
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Figure 2: Sum-rule results for the parameters — Ai (lower curves) and A 2 (upper 
curves). For each quantity, the three curves correspond to the following values 
of the continuum threshold: u> c = 1.7 GeV (dashed), 2.0 GeV (solid), 2.3 GeV 
(dash-dotted). The vertical dashed lines show the sum- rule window. 



where the first error reflects the variation with uo c and T, while the second error 
takes into account the uncertainty in the values of the vacuum condensates. When 
evolved to a high renormalization point, our result corresponds to ^{f^b) — 
(0.15 ± 0.03) GeV 2 , which is in good agreement with the value in (||) extracted 
from spectroscopy. A very similar result has been obtained by Ball and Braun 
using the same approach |24| , and by the present author using a different analysis 



based on two-point sum rules PU 



The stability of the sum rule for the kinetic-energy parameter Ai is not quite as 
good. The reason is that the condensate contribution has the opposite sign of the 
perturbative contribution. Inside the allowed parameter space for u c and T, we 
find values for — Ai ranging from -0.02 GeV 2 to +0.15 GeV 2 . The dependence 
on the Borel parameter is strongest in the region of low T values, where the 
contribution of the mixed condensate becomes very large. If we restrict ourselves 
to the region of larger T values by requiring that the condensate term be less than 
50% of the perturbative contribution, we find that the region below the hatched 
line in Fig. |2] is excluded. This leads to 

- Ai w (0.10 ±0.05 ±0.02) GeV 2 , (24) 

where the second error reflects again the dependence on the condensate parame- 
ters. It must be stressed that because of the relatively poor stability the sum-rule 
prediction for Ai is affected by systematic uncertainties that may be underesti- 
mated by the error quoted in (p4|). Keeping this reservation (which applies equally 
to previous sum- rule determinations of Ai) in mind, we note that our value in 
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Figure 3: Sum-rule results for the functions <pi{w). The width of the bands 
reflects the variation of the results with the continuum threshold (1.7 GeV < 
uj c < 2.3 GeV for the light bands, and u> c = 2 GeV for the inner, dark bands) 
and the Borel parameter (0.6 GeV < T < 1.2 GeV). 



( |24D implies an average momentum of the heavy quark inside the meson of order 
200-400 MeV, which appears to us to be a reasonable value. Clearly, our result 
is much smaller than the value — Ai = (0.52 ± 0.12) GeV 2 obtained by Ball and 
Braun indeed, we find — Ai < A2(//o) for all choices of the parameters. We 
shall comment below on the difference between their approach and ours. 

Although our main focus was to derive sum rules for the parameters Ai and 
A2, the invariant form factors <pi(w) defined in (^) may be of some interest as 
well. For instance, the combination f(w) = 3(j) 3 (w) — 2{w — l)0 2 (w) appears 
in the analysis of non-factorizable contributions to class-I non-leptonic two-body 
decays such as B° — ► D + n~ | 4"7fl . Therefore, we find it worth while to study the w 



dependence of these form factors using the sum rules in ([T5|) combined with the 
sum rule in (^). In Fig. [|, we show the results for the functions 4>i{w) obtained 
by varying the parameters u c and T in the ranges described above. For all three 
functions, we observe a mild decrease with w. We note that these results can be 
trusted only for moderate values of w, which is however sufficient for all practical 
purposes. For w 1, the non-perturbative contributions to the sum rules would 
rapidly vanish once the non-locality of the condensates was taken into account. 



3 Concluding Remarks 

We have presented a simultaneous determination of the HQET parameters Ai and 
A2 from a QCD sum-rule analysis of a 3-point correlation function containing the 
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local operator h v iTig s G^ u h v together with two heavy-light currents. Our sum 
rule for the parameter A2, which describes the chromo-magnetic interaction of 
a heavy quark, agrees with a result derived previously in Ref. pl |. However, 



our sum rule for Ai is new. It incorporates the virial theorem, which relates the 
kinetic energy of a heavy quark inside a meson to its chromo-electric interaction 
with gluons |26fl. We shall argue below that our approach is superior to that of 



Ball and Braun |24[ , who extracted Ai from the correlator of the kinetic operator 
h v (iDj_) 2 h v with two heavy-light currents. The reason is that the virial theorem 
makes explicit an "intrinsic smallness" of Ai, which is otherwise hidden by a large 
background from excited-state contributions. Our numerical results are given in 
(|23"D and (p4|) . The value of Ai implies an average residual momentum of the 
heavy quark inside a meson of order 200-400 MeV; the result for A2 translates 
into a spin splitting of m 2 B , —m 2 B = (0.60±0.12) GeV 2 , which is in good agreement 
with experiment. 

We like to add a final comment regarding the difference between our approach 
and that of Ball and Braun []24[| , who obtained a much larger value for Ai than 
our result in (|24]). The reason is that their sum rule contains a large contribution 
from a "bare" quark loop, which is 0(a° s ). Before the continuum subtraction, 
their result reads[] 

- Ai F 2 e- 2K ' T + C{T) = - | ml (qq) + 0(a s ) , (25) 

where C(T) represents the contributions of excited states, which are removed in 
the continuum subtraction. We have argued in Ref. P?J that the virial theorem, 
which relates Ai to a matrix element of the gluon field-strength tensor, does not 
allow terms not containing the gauge coupling in the sum rule for Ai. We shall 
now explain how this statement is consistent with (|25|) . 



To start with, let us stress that we do not claim that the authors of Ref. pi 
made a calculational mistake; indeed, we have checked that their result fl25"D is 
correct. What we are going to argue is that the leading perturbative term on the 
right-hand side of (|25|), together with part of the contribution from the mixed 
condensate, must not be attributed to the ground-state, but rather to excited 
states. The virial theorem helps to avoid from the start any subtleties related to 
the complicated problem of the continuum subtraction in 3-point sum rules. 

In order to explain our argument, it is necessary to go into some of the de- 
tails of QCD sum-rule calculations in the HQET. The sum rules for the matrix 
elements of local dimension-5 operators are closely related to the derivatives of 
the sum rules for some lower-dimensional operators with respect to the Borel 
parameters. Consider the 3-point sum rules for the meson decay constant (i.e. 
the sum rule for the Isgur-Wise function evaluated at zero recoil) and for the 



4 For simplicity, we do not display the terms of 0(a s ) here, since the main problem is to 
understand the origin of the leading term. 
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product £_(1)-F 2 , where £-(1) is the zero-recoil limit of a form factor defined in 



terms of the matrix element of the operator h v iTiD^h v [38|. These sum rules 



read [3/ 



F 2 e" 2 ^') +Cl (t + = A - (gg) + (t + ^ + (a s ) , 

F 2 e" 2 ^') + C 2 (t, = ^ - ^ (7f - t) + 0(«.) , (26) 

where t — 1/r and £' = 1/V are the two Borel parameters associated with the 
variables to and u/, and Ci denote the contributions to the correlators from exited 
states. As a consequence of the orthogonality of states, C\{t + t') is a function 
of the sum of the Borel parameters only. Note, however, that the second sum 
rule is not symmetric in the two Borel variables. This is a consequence of the 
fact that the operator whose matrix element defines £-(w) contains a derivative 
acting on one of the heavy-quark fields. Consequently, the function C 2 (t,t') is 
not symmetric in its arguments. Taking derivatives with respect to the Borel 
parameters t and t', we can derive a set of related sum rules containing powers 
of the parameter A. In the case of the first sum rule in (|26|), it clearly does 
not matter whether we take a derivative with respect to t or t'; however, the 
same statement is not true in the case the second sum rule. After taking the 
derivatives, we set the the Borel parameters equal, t = t' = 1/2T, in which case 
the first sum rule in (EBT) reduces to (JZUp. We obtain 



F 2 e~^ T + C x {l/T) = g - (g-g) + ^ + 0(a s ) 

A F 2 e" 2 ^ - \ C[(l/T) = + 0(0.) , 

A 2 F 2 e~ 2K ' T + \ C'{(l/T) = % + + (a s ) , 

4 An z a 



(27) 



and 



£_(!) F 2 e~ 2 ^ + C 2 (1/2T, 1/2T) = ^ - + 0(a s ) , 

loir 1 8i 

Ae_(l) F 2 e- 2 ^ T - \ d t C 2 (l/2T, 1/2T) = g - + (a s ) , 

1 ST 5 7 

A£_(l)F 2 e- 2A / r --^C 2 (l/2T,l/2T) = — + - m 2 (gg> + 0(a.) . (28) 

It is crucial that the continuum contribution C 2 {t,t') is not symmetric in t and 
t'; otherwise the last two sum rules would be inconsistent with each other. 

Using the equations of motion of the HQET, one can show that £-(1) = A/2 



[|38[| . In fact, this relation follows by comparing the first sum rule in fl28|) with 

l 
l 



the second sum rule in (p7|) , provided we identify C 2 (1/2T, 1/2T) = —\C[{1/T) 
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This relation between the continuum contributions is indeed satisfied when one 
adopts the standard continuum model, according to which 



Ci(Vr) = [1 - H"c/t)\ , c 2 (i/2T, 1/2T) = _ [1 - s 3 (u c /t)] 



(29) 

Things are more subtle for the last two sum rules in (p8|), which have the same 
ground-state contribution. Inserting the normalization condition £_(1) = A/2, 
we find that their sum agrees with the last sum rule in (2"7|). However, the only 



logical explanation for the fact that the theoretical expressions on the right-hand 
sides of these sum rules do not coincide is that the difference between these 
expressions contributes to the exited states only, but not to the ground state. 
Hence, taking the difference between the two sum rules in (|28|) leads to a sum 
rule with vanishing ground-state contribution. It reads 

ST 5 m 2 (aa\ 

{d t , - d t ) C 2 (1/2T, 1/2T) = _ - + 0(a s ) . (30) 

Let us now come back to the sum rule (|25|) for the matrix element of the 
kinetic operator derived by Ball and Braun. We can combine it with the above 
relation in such a way that the contribution of the bare quark loop, which is 
forbidden by the virial theorem, is eliminated from the right-hand side of the 
sum rule. Then the result takes the form 

- Ax F 2 e- 2R / T + C(T) = + (a s ) , (31) 

8 

where the new continuum contribution is given by C(T) = C(T) — \(d t > — 
dt) C2(1/2T, 1/2T). What we have achieved is to identify the leading perturbative 
contribution in (p5|), as well as part of the contribution of the mixed condensate, 
as a contribution to the exited states coupling to the correlation function. What 
remains is nothing but our sum rule for Ai obtained in (0), with the correct 
coefficient in front of the mixed condensate. 

Acknowledgements: It is a pleasure to thank Zoltan Ligeti and Chris Sachra- 
jda for helpful discussions. 
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